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We discuss the classical limit for the long-distance ("soft") modes of a quantum field when the 
hard modes of the field are in thermal equilibrium. We address the question of the correct semiclas- 
sical dynamics when a momentum cut-off |p| < kc <ti T is introduced. Higher order contributions 
leads to a stochastic interpretation for the efi'ective action in analogy to Quantum Brownian Motion, 
resulting in dissipation and decoherence for the evolution of the soft modes. Particular emphasis 
is put on the understanding of dissipation. Our discussion focuses mostly on scalar fields, but we 
make some remarks on the extension to gauge theories. 
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I. INTRODUCTION 

Solutions of the classical field equations in Minkowski space have been widely used in recent years to describe long- 
distance properties of quantum fields that require a nonperturbative analysis. These applications include: the diffusion 
rate of the topological charge in the electroweak gauge theory the thermalization rate of non-Abelian gauge 

fields 1^-^, as well as a wide range of cosmological problems, such as thermalization, decoherence, and structure 
formation [9|-p^. Many of these studies are concerned with dynamical properties of quantum fields near thermal 
equilibrium. 

Classical treatments of the long-distance dynamics of bosonic quantum fields at high temperature are based on the 
observation that the average thermal amplitude of low-momentum modes is large. For a weakly coupled quantum 
field the occupation number of a mode with wave vector p and frequency LUp is given by the Bose distribution 



(1) 



For temperatures T much higher than the (dynamical) mass scale m* of the quantum field, the occupation number 
becomes large and approaches the classical equipartition limit 



n{ujp) 



.0 T 



> 1. 



(2) 



The classical field equations should provide a good approximation for the dynamics of such highly occupied modes.|] 
At a closer look, however, the cogency of this heuristic argument suffers considerably. The thermodynamics of a 
classical fieldQ is only defined if an ultraviolet cut-off kc is imposed on the momentum p such as a finite lattice spacing 
a. Many, if not most, thermodynamical properties of the classical field depend strongly on the value of the cut-off 
parameter kc and diverge in the continuum limit (kc — > oo). Examples are the energy density 



e(T) 



{2Trhy 



fiLjpn{Ljp) ^ Tkl 



(3) 



and the dynamical mass gap 



^Note that for interacting quantum fields, generally, the dynamical mass gap at finite temperature is of order gT, where g is a 
dimensionless coupling constant. Hence, the validity of (^) requires weak coupling. On the other hand, rapid thermalization and 
decoherence of the soft modes requires that the coupling to the heat bath is sufficiently strong . It needs to be established 
whether both conditions can be satisfied simultaneously. 

^Indeed, the ultraviolet divergence of thermodynamic quantities was one of the motivations for Planck's conjecture in 1900 
that the excitations of field modes are quantized. 
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in the high-temperature hmit (and for fcj, ^ m*), which contrasts with the results e(T) ~ T'* and m*^ ~ in 
the full quantum theory. Remarkably, some thermodynamic quantities, such as the diffusion rate of the topological 
charge [|[|| in non-Abelian gauge theories or the maximal Lyapunov exponent Q , which is equivalent to the damping 
rate of soft thermal excitations , are found to be insensitive to the cut-off, attaining a finite value in the continuum 
limit. 

The question arises whether the value reached in the continuum limit agrees with the one obtained in the full thermal 
quantum field theory, if a complete nonperturbative calculation were feasible. Naively, one expects that any quantity 
that is independent of h in the lowest order of its coupling constant expansion — such as quantities that contain only 
the combination g^T in thermal non-Abelian gauge theories — should be reliably calculable in the classical limit. Of 
course, there will exist quantum corrections which are suppressed by powers of the coupling constant a = g^Ti, i.e. 
such a quantity will have a perturbative expansion of the form 

{g^TY{cQ + cia + C2a^ + ...) (5) 

where only the coefficient cq is calculable classically, but the coefficients Q (i > 1) can be obtained using perturbative 
techniques. 

Other questions also arise: What coupling constant g should be used in the classical calculation? How is it 
renormalized? How does mass renormalization enter into the equation for the low-momentum modes that are described 
classically? How does the classical equation "know" that it derives its justification from the classical limit of a thermal 
quantum field? How does the temperature enter into the classical equation? Which mechanism ensures that the 
attained temperature equals that of the underlying thermal field theory, if the classically described modes thermalize 
due to their interaction? 

Some of these questions were addressed in recent articles |l2| , |l8| -[20| ; others remain unanswered. Here we make an 
attempt to address the complex of questions in its entirety, applying the method of effective field theories and the real- 
time formulation of finite temperature quantum field theory. Our approach leads us to the formulation of an effective 
dynamics for the soft modes of the quantum field in terms of a stochastic, dissipative action. The stochastic forces 
describe the exchange of energy and other quantum numbers with the perturbatively treated hard modes, and the 
dissipative terms describe the eventual approach to thermal equilibrium of the soft modes. A fluctuation-dissipation 
theorem ensures that the soft modes reach the same temperature as the hard modes with which they interact. We 
perform our analysis for the simplest nontrivial example: a massive scalar field with quartic self-interaction in (3 -I- 1) 
dimensions. This allows us to evade the technical difflculties associated with gauge invariance and to make contact 
with previous work [|jl|JT|-|§. 

It is worthwhile mentioning that a similar issue arises in the Euclidean formulation of thermal non-Abelian gauge 
theories, because thermodynamic quantities, such as the pressure or screening masses, receive contributions at a 
scale (^^T)", which are essentially classical and usually cannot be calculated by means of perturbative methods. 
Recently, Braaten and Nieto (2|j2^ have shown how low-energy effective actions in Euclidean space can be utilized to 
systematically resolve this problem. The analogous problem in Minkowski space is more complicated, since dimensional 
reduction [E3j does not occur in this case because at finite temperature there is no gap in the excitation spectrum. 
In this case the effective long-distance field theory becomes essentially classical, because quantal loop corrections are 
suppressed by powers of Tikc/T <C 1 where fee is the ultraviolet cut-off of the effective theory. However, the dynamics 
retains its full (3-1-1) dimensionality. 

Our manuscript is organized as follows: In the next Section wc first review the concepts of effective actions and 
influence functionals. We then derive the influence functional for a self-coupled scalar quantum field and obtain the 
classical, stochastic equations of motion for the soft field modes. In Section HI we show in detail how dissipation 
emerges from the effective equation of motion and how the proper Markovian limit of the corresponding nonlocal 
terms can be obtained. Section IV briefly considers the generalization of our results to gauge theories and presents a 
discussion of some open issues. 

II. REAL-TIME EFFECTIVE ACTION FOR CLASSICAL SOFT MODES 

A. Preliminaries 

Our goal is the derivation of an effective action for the classical dynamics in Minkowski space of the low-momentum 
modes of a scalar quantum field $ near thermal equilibrium. Following Bodeker, McLerran, and Smilga |lq] and 
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Lombardo and Mazzitelli pC( |, we divide the Fourier components of the field $ into low- momentum modes with 
IpI < kc (henceforth called "soft" modes) and high-momentum modes with |p| > kc (called "hard" modes from now 
on) by defining^ 

$(p, t) - <i>(p, t)0{k, - IpI) + <i>(p, t)0{\p\ - fc,) = (/.(p, t) + ^(p, t). (6) 

It is clear from the introduction that the cutoff scale kc should be much smaller than the temperature T, i.e. kc <C T. 
Still kc may be chosen comparable with the dynamical mass scale m* ~ gT for the quasi-classical regime. We assume 
that all hard modes tp are, and remain, thermally occupied.^ We do not make this assumption for the soft modes 4>, 
but we assume that their occupation number is sufficiently large to warrant a classical treatment. 

We note here that this procedure is not the same as obtaining an "effective action" (in the common usage of the 
expression) for the expectation values of the quantum field in the presence of classical sources. The latter entails 
calculating the generating functional 

Z(J) ^ Jt^<^> e^(^[*l+'^« = e'^I^l (7) 

allowing for sources J of the soft modes only, and then taking the Legendrc transform 

SW 

m = W[J] - J— (8) 

with (j) = {4>) J — SW/SJ . The effective action r(0) contains all quantal and thermal corrections to the field expectation 
values. 

Our definition of an "effective action" for the soft modes i^ 

e<n[*^ = Jvip e^^[*l = Jt^^ (9) 

which excludes all quantal and thermal corrections from the soft modes themselves. The two definitions only agree 
in the limit kc 0: 

s^^^[cj,]'^'rm]. (10) 

Which approach is preferable depends on the infrared behavior of the quantum field under consideration. When it 
is well behaved, so that perturbation theory — possibly after resummation of a certain class of diagrams — provides an 
adequate description of the low-momentum sector of the field, the standard effective action T{(j>) is the best choice. 
r(0) already accounts for all thermal and quantal corrections and determines, in addition to the evolution of the field 
expectation value, all n-point functions of the field. On the other hand, if the low-momentum field modes are governed 
by truly nonperturbative physics, such as in the case of non-Abelian gauge theories at high temperature, we have to 
rely on a numerical treatment of the long-distance dynamics on the basis of 5^^^ [</)]. The major advantage then is that 
all the short-distance dynamics that can be treated perturbatively is already contained in the effective action S^g, 
thus the numerical evaluation can concentrate on the nonperturbative sector alone where quantal corrections among 
the soft modes are of minor importance while classical (thermal) effects dominate. 

The derivation of an effective equation of motion for the soft modes by elimination of the fast modes corresponds to 
a coarse graining of the quantum field. In the spirit of the general formalism of quantum statistics, the assumption of 
thermal equilibrium for the eliminated modes can be interpreted as one of minimal information about the microscopic 
quantum state of the hard modes. The general theory of coarse graining ||2^ then permits several statements about 
the nature of the effective dynamics of the soft modes: 

1. The mass parameters and the coupling constants are renormalized, and new effective interactions, which are in 
general nonlocal, will appear. 



^In the notation of ref. [^: (f> ~ ^h, ^ ~ ^s; that of ref. |^: = $<, 79 = $>. 

^This is, potentially, a weak point because the thermalization rate of hard modes is usually not significantly larger than the 
one of soft modes. By using a general form of the density matrix for the hard modes, the assumption of their thermalization 
could, in principle, be relaxed. 

^Below we will generalize the effective action to the situation appropriate for real-time evolution, which requires a doubling 
of the degrees of freedom: S^g[cj>] — » S'fp [cj>, (j>']. 
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2. The dynamics is dissipative and contains noise terms satisfying a fluctuation-dissipation theorem. The semiclas- 
sical real-time evolution is inherently stochastic, reflecting the coupling of the soft modes to the heat reservoir 
of hard thermal modes. 

3. The dissipative terms are nonlocal in time, giving rise to memory kernels stretching over the history of the 
evolution. This occurs because the heat bath needs time to respond. A change of the system (here the soft 
modes) influences the heat bath which in turn responds causally back onto the system. 

The most convenient techniques for the formulation of the problem and the derivation of the effective action are the 
closed-time-path (CTP) method of Schwinger and Keldysh pa and the influence functional (IF) method of Feynman 



and Vernon ||2^. As discussed, e.g., in ref. |11 2^ these two approaches yield identical classical equations for the 
expectation values of the coarse-grained field components (j)- Iii appendix A we summarize the basic idea of the 
influence action and its properties for a quantum mechanical system X coupled to a heat bath Q. The influence action 
S'if[0, 'P'ltf] is deflned as (see appendixes A and B) 

d^fd^.d^'J I?(^'e'(^[^l+'5'-['^^^l-s[^'l-S'-['^'^^'l)ph(</5,,<p:;i,) (11) 

J ipi J 1^'. 

where Phi^i, ^i) is the density matrix of mitial configurations of the hard field modes and 

S[^^S[<P] + S[^]+S,^,[(j,,^] (12) 

is the action for the scalar flcld. The influence action vanishes for (j) — (j)'; its variation with respect to (f)^ — {(j) — (j)') 
yields the semiclassical corrections to the equation of motion for the soft field modes (j) due to their interaction with 
the hard modes ip. 

The generating functional in the closed-time-path method is defined by introducing sources for the soft field modes 
only: 



giWcTP [ 



^■•''1 = j d^id^,d^[ /'*'2?$'e*(^[*l-^[*'I+'^'^-^'^')p(*^,$-;i^), (13) 



where here the limit ti — > — oo, t/ — *■ +oo is understood. Assuming that the initial density matrix of the full field 
factorizes into one for the hard and one for the soft modes, 

p{'^^,<^[]U) = p],{ip^,ip'-t,) ® p^{4>,,(j)'-U), (14) 

the generating functional can be expressed in terms of the influence functional as 



If the classical approximation can be made here, i.e. if loops involving soft momenta can be neglected and the density 
matrix ps{4>iT 4>i]ti) can be assumed as diagonal, the expectation values of the soft held modes are described by the 
effective closed-time path action 



CTP I 



,0']«5[0]-5[0'] + 5if[0,0';oo] (16) 



where (j) = dWcTp/dJ and </>' — dWcTp/dJ'. 

The main step to be taken in the derivation of effective classical equations of motion for the soft field modes 
are therefore: (1) calculating the infiuence action due to the hard modes, and (2) showing that loop contributions 
involving soft modes can be neglected. It is our objective in the next section to derive the infiuence action for the 
long- wavelength modes and then to discuss in the following sections the meaning of the additional terms contributing 
to the quasi-classical equation of motion for the soft modes in more detail. 



B. Influence Action 



In most cases of interest the influence action must be obtained from a suitably truncated perturbation expansion. 
This raises the question how far one has to go in the expansion to obtain all important physical contributions. Here 
we will concentrate on the case of a self-interacting scalar field with action 



4 



(17) 



where t > U describes the time of observation. Details of the derivation of Sw[4>-i(j)'] are given in appendix B. It is 
our purpose here to outline the strategy and the results obtained. As will become clear, it turns out that one must 
include all diagrams up to order g'^ in the perturbative expansion, even ii <C 1. The interaction that generates 
Sw[cl>A'\ is 



J* d^x [y^ (0-v + §'/'' 



(18) 



All non- vanishing diagrams contributing to 5'if up to second order in 5*1111 are shown in Figure 1. The diagrams 
(a) (d) are the familiar one-particle irreducible diagrams appearing in the expansion of the effective action T{(p). The 
fifth diagram (e) vanishes in the exact limit kc 0. It contributes to the coarse-grained effective action whenever 
kc 7^ 0. Other diagrams not listed vanish due to the constraint of (3-dimensional) momentum conservation. The 
inclusion of the diagrams of order is crucial, because this is the lowest order at which the real-time effective action 
develops an imaginary part which gives rise to stochastic terms as well as the corresponding real part will give rise 
to dissipative terms in the equations of motion for the soft modes. We emphasize that this property is not connected 
with the number of loops in the diagram: the two-loop diagram (c) as well as the one-loop diagram (d) and the tree 
diagram (c) contribute to noise and dissipation. 

We define the real-time propagators for the hard modes as 

Gi^Jx -y) = i{T {ip{x)^{y))) ^ G^=(x - y) 



GlUx-y) 

GtMx-y) 

GtUx-y) 



i{ipiy)ip'{x)) =G%^{x-y) 
i{ip'ixMy))^Gi^{x-y) 
i{T{^'ix)^'iy))), 



(19) 



where the superscript "fee" indicates that these propagators are defined in the restricted space of momenta |p| > kc- 
As discussed in the introduction, we assume that the hard modes are thermally populated, with a simple dispersion 
relation 



(20) 



(At high temperature where the bare mass m can be neglected, i.e. m <^ gT, choosing kc on the order of the dynamical 
mass scale gT or even smaller, one can modify this dispersion by taking into account, to lowest order, the dynamically 
generated mass (of order ~ gT) for the hard modes themselves.) 

The complex influence action Sif[(1>, (f)'] can be expanded in powers of the fields (f), (j)' . Up to order g^, this expansion 
has the following structure (see appendix B): 



Res', 



IF 



Im S] 



IF 



rt ^1 

^~ it ^^^^^^ ^ (27V)! ^^^^^ ^™ 



Rjf{x) 9{xo - yo)^e V^^ix-y) R+{y) 



(21) 



where is the amputated 2A''-point vertex function, restricted to contributions solely from the hard modes with 
|p| > kc- Here we have used the abbreviations 



R%{x) = ct>{x)'' ± cp'ix 



(22) 



Up to order g'^, the vertex functions Fj'', T'I", Tq" are given by the diagrams in Figure 1. All higher vertex functions 
vanish. Fj" is the one-particle irreducible self energy of the soft field modes as generated by the hard modes, F^" is 
the one-loop correction to the bare four-particle vertex, and Fg"" is a one-particle reducible six-particle vertex which 
vanishes in the limit kc 0. Detailed expressions for the vertex functions are given in Appendix B. 
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C. Interpretation of the imaginary part of the Influence Action 



Before one can derive classical equations of motion for the soft field modes, one has to deal with the imaginary part 
of the influence action. The general idea of its physical significance is outlined in appendix A. It has to be interpreted 
as a stochastic "external" force driving the classical degrees of freedom, i.e. the soft modes. This stochastic force is 
an important ingredient for the equilibration of the classical modes. The standard procedure [|6| is to introduce a 
stochastic variable ^ for each imaginary contribution to S'lp, making use of the identity 

e-ix"Ax = [det(27rA)]-i/2 J e-^«"^-'?+*x"«. (23) 

Here A is a symmetric matrix, x ^'^d ^ are vectors, and we have employed obvious vector notation. Applying this 
technique here we need to introduce three separate stochastic functions £,n{x), N =1,2, 3, with the stochastic weights^ 



P[S,n] = Cat exp <^ d xd y Cw(a;) 



(27V) 



(lm[r^-^])^^'^ e^(2/)|, (24) 



4iV2 

where Cn are appropriate normalization constants. The influence functional then takes the form 



J V^,V^2V^3 P[^i]P[C2]P[^3] e'^-[*'*'^«-l (25) 



with the new stochastic influence action 

^iF a, 6, 6] =Re (5iF [(/.,</.'])+ f d^x Rj,{x)iN{x). (26) 

N=l 

This expression forms the basis for the classical equations of motion for the soft modes 0(x). 

D. Classical Equations of Motion 

As we argued in the introduction, the dynamics of the soft modes under the stochastic action 

5[</)]-5[0']+5if[^,0';^] (27) 

is governed by (quasi-)classical physics, if we assume that the soft modes are sufficiently populated. As we will 
argue later this assumption will be justified by the full dynamics which we will derive and discuss in the following. 
In addition, as sketched briefly in appendix A, we note that possible initial quantum mechanical correlations in 
the system will decohere due to the stochastic forces acting randomly on the soft modes and so destroying those 
correlations (or phases) with time. We will not address this issue here, as in the present case we are convinced 
that the soft modes behave already classically, and thus refer the reader to the extensive literature on decoherence 
[ p^ , pO 29-3^. Loop contributions to the effective action rcTp[^] from the soft modes (under the influence of Syp) 
are not necessarily negligible, but at high temperature they are dominated by thermal fluctuations, not quantal 
ones, which are relatively suppressed by a factor tikc/T . Thermal loops correspond to tree diagrams with thermally 
populated external legs. They can, therefore, be accounted for by averaging the classical equations of motion at the 
tree level over a thermal ensemble of initial states. 

The classical equations of motion for the soft modes are obtained by a generalized stationary phase approximation 
(see appendix A) for the effective action defined on the time contour path C by first varying ( p7| ) with respect to 
(/)A = (0 — (/)') and then setting (j) = (j)' . 

Making use of the relations 



^We have chosen the index of the stochastic field to match the power of the field operators in the stochastic action. In the 
notation of ref. ^2 £,1, £,\ C2; in that of ref. Q. ^1 ^ —rj, ^2 ^ — 
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we find: 



5R%{x) 



N^{x)^-H{x~y), 



0=0' 



= 0, R%(x)\ =24>{x)'', 



(28) 



52 

mMx) + rri^cj)(x) + —4>^{x) + \^ 
6 ^-^ 



\N-1 



N-- 



'^{2N~l)\ 



d'^yRe 



(29) 



Ar=l 



We first note that tlie resulting equation of motion for the soft fields is explicitly causal as the real part of the 
influence action (|2^) contains a step function assuring causality, i.e. 



0{x„ - yo)Re T'^n{x - y) 



(30) 



defines just the response function (to lowest non-vanishing order) of the hard modes on the soft modes and thus 
corresponds to the familiar Kubo relation of systems near equilibrium generalized to the present case. 

Secondly, the real part of the infiuence action is divergent and must be renormalized. Since the ultraviolet diver- 
gences of the vertex function are independent of the low-momentum cut-off and of the temperature, they 
can be absorbed in the standard counter terms for mass and coupling constant renormalization. We can, therefore, 
replace the divergent T^'j^{x — y) by the regularized expression 



t^J^ix -y)^ r'^J^ix ~y)- G2n5{x - y) 



(31) 



where G2N subtracts the temperature independent short-distance singularity of V^"^ 



In our specific case here, 



and r^'' require renormalization. Tq" does not because it is finite and even vanishes in the limit fc^ 0. Hence, 
TO and g denote the renormalized mass and coupling constant, respectively, and f jat the renormalized vertex 
functions. The tadpole contributions of Figure 1(a) and 1(b) to r2 contribute thermal terms to the effective mass, 
changing m to m* . The thermal effective mass to* will also depend on the momentum cutoff fee as only the mean 
field contributions of the hard modes are accounted for. Equation ( p9| ) is graphically represented in Figure 2. The 
terms on the right-hand side represent stochastic source, mass, etc., terms for the soft field modes. In particular, the 
term £,i{x) indicates that the soft modes do not obey a source- free equation, rather, they are in constant contact with 
randomly fiuctuating sources generated by the hard modes that have been integrated out. 

The equation of motion ( p9| ) resembles a Langevin equation analogous to the situation of quantum Brownian 
motion. However, as it stands, it is clearly nonlocal in time, and one has to identify the appropriate Markovian limit 
(or approximation) of this equation. Indeed, the question arises whether a Markovian limit exists at all. The answer 
to this question will have a bearing on the detailed understanding of dissipation in quantum field theories. Clearly, if 
a Markovian limit exists, it will provide a much simpler equation which can be applied for practical purposes. 



III. DISSIPATION 
A. Origin of Dissipation 

In this section we elucidate the dissipative character of the resulting contributions to the semiclassical equation of 
motion (|2|) deriving from the real parts of the influence action at order g'^. Only the contributions of Re ^S'jp''^ , 

Re {s[^^ and Re 

(•^IF^) will be responsible for dissipation, whereas the other contributions only renormalize the 

effective mass parameter up to order g'^. With this in mind we concentrate on their influence on the equation of 
motion ( p9|) of the soft modes (f>{x): 
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\i—c,d,e 



6(1) Aix) 



(32) 



According to the general fluctuation-dissipation theorem, the three terms on the left-hand side lead to the dissipation 
of energy stored in the soft modes (the "system") to the hard modes (the "bath") balancing the energy gained 
continuously from the fluctuating force terms on the right-hand side, which were obtained from the imaginary parts 



of the influence functional Im ( 5, 



Ac) 
'if 



, and Im (s*};') . 



From Appendix (B) we immediately find for the three dissipative terms in the sum of eq. (p2[): 



dRe 5, 



64 




6Rc 




6(f} 


a{x) 


6Re 





6(1} a{x) 



i=0 



4a=0 



(G5.^(i,T)3-G^^(i,r)3 
d'iC-g') (G^^(x,r)2-G^i5c,r)2) 



X — X, i — r) 



x,i 



(33) 



x,t-r)^ 



where we have set xq = t. 

Obviously, the contributions of the hard modes have a rather simple structure: The integration kernels are causally 
weighting the difference between powers of the real-time propagators G^" and G< . For convenience we define the 
memory kernels Ai^^^: 



1-4 

6^ 



G^=(x,r)^-G^=(x,r) 



M^'\i,r) = (G^=(i,r)2 - G^< (A,r)^) 



(34) 



X(")(x,t) 
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G^=(x,r)-G^^(x,r; 



To further evaluate the dissipative terms ( p3| ) one has to insert an approximation for the time dependence of the soft 
modes (/)(x — x, < — r) under the integrals. Previous authors [p[ p^Jl9| , |28[ | have expanded around t up to the first 
order gradient correction, i.e.: 



x,t-r)^ 



(x,i)" - nT(/)(x,t)"-V(x,t) 



(35) 



This amounts to a quasi instantaneous approximation of the soft modes in the Heisenberg picture. If this approxima- 
tion is made, the dissipative terms ( ^3| ) vanish, unless an explicit width T is introduced]^ in the propagators G^'^ ^ of 
the hard modes §,0. 

We believe that this procedure obscures the true dissipative character of the memory terms. In fact, as we will 
see below, the Fourier frequencies of the soft modes are essential parts of the correct strategy for revealing the close 
connection of the expressions (^3|) with the more familiar on-shell scattering rates of soft modes on hard particles. 
Only if the soft modes oscillate periodically with some given frequency, can they scatter on the hard momentum modes 
or, if the frequencies allow for this, even produce hard particles. In a certain sense, the soft modes need time to scatter 



^The inconsistency of this approach becomes apparent when one compares the memory time of the dissipative kernels with 
the characteristic oscillation period of the soft modes {{m*^ + k^)~^^^). From (^) the extension in past time of the kernels 
can roughly be estimated to be of the same order, i.e. k^^, as it sets the the typical scale. (The memory time in a scattering 
process is given by the uncertainty principle as Tmem ~ 1/ ( A_B) , where {AE} denotes the averaged energy transferred in the 
reaction. For the soft modes the energy transferred cannot be larger than kc.) Hence the soft modes will oscillate entirely over 
the extension of the memory kernels. 
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and again come on shell. The instantaneous approximation ( |35| ) clearly does not account for the oscillatory behavior 
of the soft modes and hence cannot properly establish the correct Markovian limit of the expressions (|33| ) . The main 
objective of the remainder of this section is the derivation of a well-defined Markovian form of the dissipative terms 
( P3| ) that takes the oscillatory character of the soft modes in account. 

The soft modes will, of course, oscillate with a certain spectrum of frequencies. This is true also for the spatial 
Fourier modes 0(k, t). In a linear approximation one expects, however, that they oscillate with some pronounced 
frequency 



Ek ~ Vm*2 + k2 



(36) 



This should be rigorously valid in the weak coupling regime when the soft modes propagate almost freely. Without 
further specifying the single-particle energy Ek , we can approximately express the soft modes at earlier times {t — r) 
by their value at the time t in the linear harmonic approximation 



(j)(k, t — t) K 0(k, t) cos EkT 



(j)(k,t)-^ sin EkT 
Ek 



(37) 



Note that the 0(k, t) contribution turns out to be essential and cannot be neglected. Comparing (^7|) with ( pSj ) one 
immediately recognizes that the harmonic approximation (^7|) contains much more information about the past than 
the instantaneous approximation ( p5[ ) . Equation ( |37| ) can be considered as a Markov approximation to the soft field 
modes in the interaction picture, where the free oscillation frequency has been extracted. 

To reveal the full meaning of the three dissipative terms ( |3^ ) in the equation of motion (^) of the soft modes 
we insert the harmonic approximation ( ^7| ) and further evaluate the three terms approximately. Our strategy is 
completely equivalent to the quasi-particle approximation used in the more familiar kinetic theories when evaluating 
the collision term in transport processes | |3^ . A certain mode (quasi-particle) is assumed to propagate with a specific 
frequency, the quasi-particle energy. This approximation is valid if the spectral function of the quasi-particle is sharply 
peaked at this frequency. Moreover, pursuing the analogy with standard kinetic theory, the equation of motion ( |3^ ) 
becomes instantaneous as all soft quasi-particle modes (/)(k, t) are evaluated at the same time. We thus obtain the 
Markovian limit of eq. (^) when evaluating the time integral over the memory kernels defined in ( p^ ) with the soft 
modes propagating according to (^). 

We write (IS^) as 



SRe [S^^) 



S4 




SRe 




54. 


A (a;) 


5Re 


(40 



(Pk 

(2^ 



,=0 



^T(")(k,t) 



^■''T('')(k,t) , 



i=0 



4a =0 



(27r)3 



(38) 



with the notations 

T(")(k,i) =6'(fcc- |k|) / dr0(k,i-T) Al(")(k,T) , 



T<-^Hk,t)= I -^eik,-\k-ki\) dr 0(k - ki, f - r)0(ki, i - r) M^'') (k, r) , 



(39) 



T(^)(k,t) 



/fee j3jL j3i, i-oo 
^^^0(fc,-|k-ki-k2|) J dr0(k-ki-k2,t-r)0(ki,t-r)0(k2,i-r)X(^)(k,r) 



where the momentum integrations for all soft modes are restricted by the cut-off momentum kc. We can now insert 
the approximate form ( |37| ) for the soft modes. It is convenient to first perform a Fourier transformation of the memory 
kernels 



(40) 
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and insert the expressions for the hard propagators of Appendix B. One finds 



X { [(1 + ni)(l + n2)(l + 713) - nin2n^] Si^J - - W2 - 0^3) 

+ [(1 + rti)n2(l + ns) - 7ii(l + 712)713] 6{uj - lui + uj2 - UJ3) 

+ [77,1(1 + 7l2)(l + ns) - (1 + «l)772«3] S{UJ + LUl- UJ2- ^^s) 

+ [71l?l2(l + 7^3) - (1 + + "-2)7^3] S{lU + UJ1+UJ2- t^s) 

+ [(1 + ni){l + 772)713 - ril7l2(l + 773)] 5{lU - LUl- UJ2+ UJ3) 

+ [(1 + 71l)r72773 - 7T,l(l + r72)(l + 773)] S{lU - LUl + UJ2 + UJ3) 

+ [ni{l + 772)773 - (1 + 771)772(1 + 773)] 5{uj + wi - W2 + UJ3) 

+ [771772713 - (1 + 77l)(l + 772)(1 + 773)] (5(ci; + CJl + ^2 + W3)} (41) 



where uji = 01^-^,102 = ujq2,^3 — t^k-qi-q2 E^nd 77^ = 77(wq^), and 

M^'\Kcu) = -Is' (^)/^^ _^li_0(|k-q|-fc,) 



(27r)3 OJl LU2 
X { [(1 + 77l)(l + 772) - 771772] S{uj - UJl - UJ2) 
+ [771(1 + 712) - (1 + 711)772] (5(tJ + Wl - W2) 
+ [(1 + 771)772 - 77l(l + 772)] (5(tJ - CJl + CJ2) 

+ [ni772 - (l+77l)(l+ 772)] (5(w+cJi+a;2)} (42) 



where uji — LUq, 102 — "-^k-q, and 



A^(^)(k,..) = -^.94(7)-0(|k|-fc,). 

{[(1 + n(uJ^,)) - 77(a)k)] 6(uj - uj^) + \n{uj^) - (1 + n{uJ^,))\ 6{uj + uj^)} ■ (43) 

The memory kernels A^'*-'(k, w) have a very simple structure. They all are purely imaginary and satisfy the relation 

>l'')(k,w) = -A^W(k,-w) . (44) 

Moreover, the kernels are directly related to the discontinuities of the respective self-energy insertions of the hard 
modes. For a detailed discussion and interpretation of such discontinuities in finite temperature field theory we refer 
the reader to the article by Weldon [Q. From the expressions (^T[-^3|) we see that each kernel is composed of two 
contributions: the direct part r£;(k, w), or the "loss" term, and the inverse part ri(k, w), or "gain" term. Following 
iQ , we write 

iMik, Lu) (xTdik, uj) ^r,{k, Lo) , (45) 

where F^ and F.^ , respectively, are positive and can be identified with the decay rate of a quasi-particle with momentum 
k and energy uj into the allowed open channels or the production rate due to the inverse processes. In addition, the 
detailed balance relation 

^^=exp(c./r) (46) 
F.,(w) 

follows immediately if the distribution function 77(0;) of the hard modes is assumed to be in thermal equilibrium. 

This identification represents an important step toward understanding the dissipative nature of the terms (|3^). 
Intuitively, we expect dissipation to be closely related to the on-shell scattering rates of the soft modes when interacting 
with the hard thermal particles. We have to manipulate the terms further in order to fully expose this connection. 

We now insert the approximation (^^ for the soft quasi-particle modes into the semi- infinite integrals T^'^\\s.,t) 
given by (^). In order to evaluate the time integral we express the memory kernels by their Fourier transform 
A^*^*^(k, w) and make use of the relation 

' 1 1 

dr e^^T-'^)^ = i lim ^ = iP -|- nSh - uj) . (47) 

e^0+ (7 — Cj) + 76 J — UJ 
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The remainder of the calculation is straightforward but cumbersome. One has to express the sine and cosine contribu- 
tions by exponentials and collect all terms. Relation (^4|) helps in identifying certain pairwise identical contributions. 
The final result reads 

T(^)(k,t) = 0(k,i)0(fc,-|k|) / ^^M^'\k,iu)P-^ + <P{k,t)9{K~\k\)-^^M^'\KEk) , (48) 

J 27r Lk - Lu 2bjk 



TW(k,t) = y'''|^0(fc,-|k-ki|)|0(k-ki,t)0(ki,i)i j '^^M^^\k,uo) 



xiP i + P- ^ 



>(k-ki,t)0(ki,t)^^-l-^ j ^zA^W(k,..) 

x[-p^ — — + p- ^ 



■ r ~ ki, i)0(ki, f)^^^0(fc, - ik - kii 



X I IJ 



iA^('*)(k,Sk-ki +SfeJ+iA^^'''(k,£;k-ki -Sfej) , (49) 
T(^)(k,t) = /'^ "^'(2'^^' ^(fcc - |k3|) |0(ki,O0(k2,^)0(k3,t) \ y'^(z)A^(^)(k,c.) 



X ( P— — ^— h P— ^ ^ ^ h 2P- 



-E'l + El + £^3 — w -El + E2 — Ej, — Lu El — E2 + £^3 — w 

<^(ki, i)0(k2, t)^(k3, t) 1 ^ * A^(^)(k, C.) 

-i"^; ^ ^ ^ i^- — ^— + 2P- ^ 



-|- E2 ~\~ — LU E\ -\- E2 — E^ — cu El — E2 ~\~ E'^ — uj 

0(fc, - |k3 1 ) 1 0(ki , t)0(k2 , t)0(k3 , t) — 

M'^''\k, Ei+E2 + E3)+ iM^^^k, ~Ei-E2 + E3) + 2iM^^\k, Ei ~ E2 + E3) 

+<^(ki, i)0(k2, t)0(k3, t) ^TT^Tir 

X (iM'-"^ (k, -El - £2 - ^3) + iM^"^ (k, E1+E2- E3) 

+2tM'^'''>{k,Ei~ E2 + E3))Y (50) 

where k3 = k — ki — k2 . 

One notices that the contributions with an odd power in (p contain the on-shell reaction rates ( ^ifji^ ). Thus we 
expect them to describe the dissipative aspects of the equation of motion for the soft modes as energy is transferred 
by scattering (or production) processes between the soft and hard modes. Since this interpretation is most obvious 
for the first term T^'^\ let us start by discussing this term in more detail. Fourier transforming the equation of motion 
( ^ ) for the mode k, we find from (38) that Re(S'|p'') contributes on the right hand side as 

4>{k, t) + k^C^{k, t) + TO20(k, + y y" ' ^^^^(fcc - |k - ki - k2|)</>(ki, i)0(k2, t)</)(k - ki - k2, t) 

f duj iM^''\k,uj) ■„ ^iM^''\k,Ek) , 

+HK t)pj - + ^(k, t) + . . . = 6(k, t) + . . . . (51) 

The first additional term can be readily interpreted as momentum dependent contribution to the mass term (after 
proper renormalization) . The second term has the familiar structure of a velocity dependent damping term, i.e. 
ri4>{k,t). From (^) and the detailed balance relation ( ^ ) we find that the friction coefficient 

!ii;f|^>0 (52, 
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is in fact positive. It is exactly related to the rate (r^ — F^) with which a small disturbance relaxes towards equilibrium 
by coupling to the various open channels . The damping rate for the soft modes in this particular channel follows 
immediately as 

7W(k) = i,7(^)(k) , (53) 
as the left hand side of eq. (^l|) has the well known solution of a weakly damped oscillator 

<^(k,t) ^ e-''<^'('')*/2cos(i?2_i^W(k)2)i/2i . 

Generalizing this observation we expect a similar interpretation of the two other dissipative terms, T^''-' and T*^'^-', 
to hold: (1) The principal value contributions give rise to a (generally momentum dependent) vertex renormalization 
and to a new (again momentum dependent) six-point vertex; (2) the on-shell contributions can be associated with 
damping processes involving two hard particles and two soft modes or involving one hard particle interacting with 
three soft modes, respectively. In Appendix C we present a further analysis of the contributions from Re(S'|p'') and 
Re(S'jp^), where a bridge is cast to the evaluation of the discontinuities by means of soft quasi-particle propagators 
rather than soft fields. In particular the amplitude squared is easily identified with the quasi-particle density of the 
soft modes, i.e. 

|0(k, t)f = 0(k, t)0(-k, i) « N{k, i) « -L (1 + N{k, tj) (54) 
in the quasi-classical regime where the soft modes are strongly populated. One finds: 

A<'''>(k + ki,£t + £t,) + iA<''''(k + k,,Ei.-£i,)) , (55) 



X 



/^)(k,t)«-|^ / i,^i0(ki,t)n0(k2,t)p 



3 /■'^'= (fiki(Pk2 

AE^J (2^)6 

X (i7W(")(k-fki+k2,Sfe-f Sfe, -f SfeJ 

+2iM^'-\k + ki+k2,Ek+Ek,-Ek,)) . (56) 



These are the obvious generalizations of (^2|) and correspond to the friction coefficients for the possible on-shcU 
scattering or production processes among two soft and two hard modes and three soft modes and one hard mode, 
respectively. We note, however, that these coefficients do not need to be positive for all times. From ( ^ ) and (^) it 
follows that Ai is positive for a; > 0, whereas M = ~Ai{\uj\) becomes negative for uj < 0. Hence the second integrand 
of (^5|) as well as the last two of ( |56| ) can become negative, so that in principle "anti-damping" can occur in some 
cases. This is not a conceptual difhculty. Suppose one would prepare an initial system where a few modes are much 
less populated than all the others. These special modes should grow in magnitude by scattering processes either due 
to the direct (/I'^-term in (^2|), or by scattering processes including other soft modes and hard particles. Energy is then 
transferred from the other soft modes to those few and into the reservoir of hard particles. 

Summarizing this section, we succeeded in properly identifying the origin of dissipation from on-shell scattering 
processes of soft quasi-particles. Dissipation corresponds directly to the imaginary parts of the self energy insertions 

ReF2^ (k, u) in the response term of eq. ( p9| ) — which should not be mixed up with the imaginary parts of ( [2l| ) 

of the effective action — as one would have expected from standard finite temperature field theory. Furthermore 
instantaneous, but momentum dependent dissipation coefficients have been derived in the proper Markovian limit of 
the effective terms (E^). 



B. The limit 



To draw comparison with other work 12 1^ it is useful to consider the limit kc for the scalar and massive 
theory: 
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0(x,t) — > 0(i); 0(k,t) (bimnyS^ik) . (57) 

The analysis follows in complete analogy to the preceding section so that we skip the intermediate steps and state 
the final results {kc 0): 



V 



^:^^M^'H0,Eo)>0, (58) 



2E, 



= Htf ^ I yfrj A (1 + 2"K)) ^(^0 -oj,) > , (59) 
AZhQ J^^^Q [2tt)-^ ui^ 

f]^''\t) — ^0 . (60) 

Here Eq denotes the infrared limit of the dispersion relation E'k of the soft modes. The first dissipation coefficient 
due to the "sunset" graph (Figure Ic) contains contributions due to scattering as well as particle production. The 
second coefficient is only nonvanishing when producing hard particles, i.e. when the low momentum mode (j)(t) 
possesses Fourier frequencies that allow for such a dissipative channel. The last term vanishes as the hard propagator 
is restricted to momenta |q| > kc = 0. 

The dynamical mass (m*)^ to order is given by the tadpole diagram (a) in Figure 1: (m*)^ = g^T^/2A, which 
is the dominant contribution to the mass in the high temperature limit, i.e. when m <C m*. Inserting this for the 
dispersion relation of the energy, i.e. Eo = m*, the damping coefficient 77^^) is exactly related to twice the damping 
rate of plasmons in theory [p5|-^ : 

= 27(m*,0) = -f^ , (61) 

dzv 24:TT 



as one would have expected by our arguments given at the end of the last section. 

We admit that our approximate treatment of a single harmonic approximation ( ^7|) for the soft mode may fail when 
(j){t) contains a whole set of frequencies, e.g. if one thinks of unstable phase transitions where (/)(t) acts as the relevant 
order parameter. In this case, however, one can work with the explicit temporally nonlocal effective terms which 
one can easily write down from the expressions of the memory kernels at k = 0. Such a procedure was undertaken 
by Boyanovsky et al. ]l^ . At zero temperature the authors examined in detail the damping associated with the 
contribution of the 'fish-graph' (d) in a more self-consistent treatment by resumming its influence on the evolution 
of the (ensemble) averaged order parameter {(j){t)). They found that the dissipation was strongest in the Goldstone 
sector when (nearly) massless particles can be produced. Such a reasoning follows also immediately from our result 
(p9|) as then rj^'^^ becomes nonvanishing even for kc 0. As outlined above, we share the authors' point of view 



that an instantaneous approximation as assumed in (35) is not valid and cannot correctly account for the dissipation 
contained in the nonlocal expressions, such as (p3|). 

C. A remark on the fluctuation-dissipation relation 

The dissipative terms and the noise terms are related by fluctuation-dissipation relations. From Appendices A 
and B we see that the memory kernels A4 = Dji (the response functions) and the noise kernels v = Dj are closely 
related (cf. their definitions (A21) or (Bll, B12, B13| )). On general grounds intimate relations among these kernels 



exist, known as the generalized fiuctuation-dissipation relations (FDR). However, their explicit realizations are quite 
complicated. For a detailed derivation which can be readily applied to our case we refer to [|lO|. Because of the 
symmetry relations A1(k, uj) = — A^(k, — cj) and i/(k, lu) — i^(k, ~lu) one defines a new kernel 7 as 

M{t) = ^7(r) . (62) 
Then a relation among 7 and u can be obtained in the form 

i/(t) = / dsK{T- s)'y{s) , (63) 



which is just the generalized FDR. In order to see how K can be derived, we make use of (^6|) together with ( |62| ) to 
write 
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7(c^) - -M{uj) cx - (rrf(k, Lj) - T,{k, uj)) . (64) 

UJ iO 

For the noise kernels a similar identification holds (compare the explicit forms given in Appendix B) 

v{u) ^TdO^,uJ)'¥V^O^,uJ) . (65) 
(The Fourier transform of the noise kernel is thus always non- negative.) Because of (^^, one finds that 

T.-( \ iM 1 cxp(ntj/T) - 1 

^("^^K^ = ^ exp(MT) + l - ^^^^ 

Now the important point to stress is that the frequencies involved for the soft modes are much smaller than the 
temperature T, i.e. <C T for kc <C T . Hence only those low frequencies are relevant for the kernels AA, v and 7, 
as also only those noises containing small frequencies can influence the dynamics of the soft modes. For v = I this 



follows in particular if one rewrites the Gaussian distribution of the noise, e.g. (A29), in the frequency domain as 

PUu)] - 1 exp {-\ J dLoaio)rH-co)a-co)^ . (67) 

If one thus considers only the relevant low frequency part of the kernels, K becomes effectively constant and thus 
local in time: 

K(s)'^2r/n, (68) 

the famous Kubo identity. Such a simple relation always occurs if the system behaves classically |^^. In fact, as we 
have argued before, the soft modes should satisfy this condition. 

The FDR ensures that the soft modes approach thermal equilibrium precisely at the temperature T of the hard 
modes, when they evolve under the equation of motion (|3^). The presence of these terms is essential, as it forces the 
soft modes to thermalize at T independent of their initial configuration. In particular, the noise terms will continuously 
'heat' the system (for field theories see rcf. |^l|]) whereas the dissipative part of the response function counteracts. 
Equilibrium is achieved when the system has thermalized to the temperature dictated by the bath. This aspect is 
not important if one only wants to study the time evolution of thermal configurations of the soft modes. However, 
for the study of configurations far off thermal equilibrium the presence of both noise and friction terms is crucial, as 
it ensures that the soft modes will become populated thermally, i.e. their amplitudes become large. This justifies our 
basic assumption. 

This remark bears relevance to the issue whether thermal masses can suppress chaos in certain nonlinear classical 
field theories as suggested by Blaizot and lancu |^]. They argued that small amplitude, homogeneous field configu- 
rations in non-Abelian gauge theories oscillate regularly due to their thermal mass, and only large amplitude waves 
become chaotic. Our above result shows that soft modes do not remain of small amplitude when stochastic forces are 
included in their description, rather their amplitude always approaches the thermal limit. It would be interesting to 
study this phenomenon, e.g. in the case of two massless scalar fields (j) and Xi which are coupled through an interaction 
of the form g^^^x^- This theory exhibits chaos in the naive classical limit. Whether the effective stochastic classical 
theory for soft field modes remains chaotic in spite of thermal mass generation remains to be seen. 



D. Memory effects 

As discussed at length above, the real part of the influence action yields nonlocal, dissipative terms in the effective 
classical equation of motion (|3^). For further exploration of the nature of dissipation we made the Markov ansatz 
to calculate the dissipation coefficients. If the strength of the soft modes, \/|0(k, t)p, also changes rapidly with 
time (e.g. in the strong coupling regime) the Markov limit may not be valid. This can happen as the memory 
kernels M{t) have support in the past roughly given by the scattering time of the individual processes. To avoid such 
interference effects one has to choose the appropriate cutoff kc large enough to restrict the temporal extent of the 
kernels. This becomes relevant in the case of massless theories in the strong coupling regime. However, the validity 
of the quasi-classical behavior of the soft modes is then no longer warranted. 
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E. Infrared behavior 



We restate the quasi-classical equation ol motion for the soft modes </)(k, t) in their approximate instantaneous 
form: 

9^ _L „W 




® + E « E ® ' (69) 

i=c,d,e N=l 

where (E) denotes a convolution in momentum space. The coefficients fii generahy depend on k and t, as well as on the 
momentum cut-off kc- Explicit calculation shows that the leading cut-off dependence of /ij derives from the thermal 
one- loop contribution (see also [39|]): 

x2rr2 ~2i 



^^ik — ' for m = . (70) 



The vacuum contribution vanishes quadratically for small k^. The same is true for the one- loop correction to the 
couphng constant 

^ (rf) k;^o ZcfT _ 2K_\ rh^ 3^/6.3 , ik,\ 9g^k, 
where m*^ — fh? + Mi(0) is the effective thermal mass. As discussed in |18 , the cut-off dependence of /^i is exactly 



' " loTrm* V 7rm* / Stt V ttJ / tt^J 



balanced by the dynamically generated mass due to the self-interaction among the soft modes, if these are also in 
thermal equilibrium. The same holds for the other fcc-dependent constants in (|69|). If T^Yi^ contains an infrared 
divergence, as often happens in massless quantum field theories, it is necessary to "resum" the thermal self-energy by 
including it into the propagator of the hard modes [p5|,p8[. (See also Il2j for an extended discussion of this procedure.) 



IV. CONCLUSIONS AND OUTLOOK 



In summary, we have succeeded in deriving a consistent set of temporally local transport equations for the long- 
distance modes of a self-coupled scalar field. Our approach is based on the analysis of time scales showing that the 
soft field modes, at high temperature and weak coupling, oscillate on a time scale shorter than the characteristic 
damping time. A temporally local (Markovian) equation can therefore only be derived for the Fourier components of 
the soft field modes, i.e. the occupation amplitudes, but not for the field itself. In the spirit of quantum field theory, 
this corresponds to the use of the interaction picture rather than the Heisenberg picture. 

As a result of our modified approach, we obtain finite values for the dissipation coefficients without the need to 
introduce an explicit damping for the hard field modes. The dominant damping mechanism for the infrared field 
mode (fee = 0) is generally given by the two-loop "sunset" diagram, corresponding to emission or absorption of the 
soft mode on hard thermal quanta. This differs from the analysis of Gleiser and Ramos who found that soft-mode 
damping is dominated by scattering on hard particles in the traditional approach of Morikawa Q . 

How do our results generalize to other field theories? Most interesting quantum field theories, especially gauge 
theories, contain three-particle vertices rather than, or in addition to, the four-point coupling investigated here. In 
these cases, there exists the standard one-loop diagram, shown in Figure 3a for the (/)'^-theory. The real part associated 
to this diagram has an imaginary part (in the (A;, a;)-representation) in the space-like domain, corresponding to the 
fluctuations in the mean field at space-like distances induced by hard thermal particles. While these fluctuations do 
not directly produce damping of on-shell soft fleld modes, they provide a mechanism for such damping, when coupled 
back into the nonlinear equation for these modes. Graphically, to lowest order, this mechanism will generate a higher 
order diagram as shown in Figure 3b. This corresponds, similar to our arguments presented in Appendix C, to the 
imaginary part of the two loop diagram (also shown in Figure 3b) when cut through the hard thermal loop and the 
soft fleld line. This analysis of mean fleld damping in the case of three-point vertices is completely analogous to that 
performed by Braaten and Pisarski |5^] for the damping of soft modes in non-Abelian gauge theories. 

Non-Abelian gauge flelds at finite temperature would present an interesting application of the methods discussed so 
far, because their infrared sector remains non-perturbative even at high temperature due to the absence of perturbative 
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screening of magnetostatic fields. This phenomenon has been widely studied in the Euclidean formalism, and its 
resolution by a sequence of effective actions has been proposed [^ , p2| . As mentioned in the Introduction, the 
situation is more delicate in the real-time formalism because of the lack of a gap in the excitation spectrum. 

The emergence of a finite correlation length in the magnetostatic sector is essentially a classical phenomenon which 
occurs at the momentum scale g^T p^-^. It should be fully contained in the classical effective action obtained by 
integrating out hard modes with momenta |p| > kc > g^T. In fact, our discussion in Section 11. E suggests that we 
should choose kc ^ gT in non-Abclian gaug e theories because the thermal gauge boson mass to* is of order gT. This 
mass itself is of semiclassical origin [[l7|-[l9| , but requires the use of the Bose distribution for the occupation number 
of hard modes and hence is sensitive to quantum physics. 

The problem in gauge theories is that one cannot employ a low-momentum cut-off kc because it is not gauge 
invariant. A lattice formulation provides a gauge invariant cut-off through the lattice spacing a, roughly corresponding 
to kc — but it violates rotational invariance. As pointed out in ref. [ p^ this necessitates the introduction of 
rotationally non-invariant counterterms. An infrared cut-off scheme that violates neither gauge nor Lorentz invariance 
is therefore preferable. Such a scheme, based on the proper time representation of the one-loop effective action, has 
recently been proposed ||5C| ]. 

We found it convenient to formulate the effective equations of motion ( |5l| ) for the soft field modes in momentum 
space. The question arises whether these equations take a local form when written in coordinate space. This issue 
becomes important for gauge theories, where the lattice regulated coordinate space representation appears most 
natural. In order to answer this question, it is necessary to consider the momentum dependence of the kernels 
A^*^*^ (fc, w) in (51). In the extreme infrared limit kc ^ Q only the term A^^'^) survives. The momentum dependence 
of this damping coefficient r]^'^'^ has been analyzed by Jeon and by Wang and Heinz , who found that it varies 
only weakly over the range < |k| < m* . Outside that range, however, a strong momentum dependence l/|k| 
is seen. A complete analysis of the two other damping coefficients rj^'^'^ and r]^'^^ has not yet been performed, but 
a strong momentum dependence appears likely. On the other hand, these terms are subdominant in the infrared 
domain. Given that Ek is also approaching a constant for |k| ^ m*, it may be possible to obtain a quasi-local form 
of the equation for the soft modes at distances larger than 1/to*. Clearly, a more complete analysis is required, 
also with regards to the important issue whether the noise terms can be approximated as local white noise, before a 
detailed understanding of thcrmalization of soft modes in quantum field theory is achieved. We hope to address these 
questions, as well as the issue of the numerical implementation of the equations derived here, in our future work. 

NOTE ADDED 

When finishing the manuscript the authors became aware of a recent manuscript by Boyanovsky, Lawrie and Lee 
on "Relaxation and Kinetics in Scalar Field Theories" js^ , where a linearized equation of motion for the ensemble 
averaged field ($(x, t)) is derived containing the contribution (a) and (c) of the real part of the effective action (cf. 
Fig. 1). Their stated damping rate coincides with our result (|6l|). 
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APPENDIX A: SHORT REVIEW OF THE INFLUENCE FUNCTIONAL TECHNIQUE (FEYNMAN) 



A quantum mechanical system X (described by the variable x) interacts with a bath Q (described by the variable 
q). At some specified initial time t — ti the combined system (X U Q) is described by the full density matrix 

Pxuq{U) ^ pxuQ{xz,qz,x'^,q'^;U). (Al) 

In the Schrodinger representation the (full) density matrix evolves in time according to 

PxuQ{x,q;x',q';t) 

dxidqidx[dq[ {U{x,q;Xi,qi;t, ti)pxuci{xi,qux[, q[; ti)U'^ {x[, q[; x' , q' -,1,^)} , (A2) 



where the evolution operator reads in the path integral representation 



U{x,q-x,,q,;t,U)^ I Vx T Vq e^^'^^'^^ (A3) 



with 



S[x, 9] ^ ds C{x{s), q{s)). (A4) 

By formally integrating out the bath degrees of freedom q one obtains an effective interaction Scff[x{s), x'{s)] describing 
the evolution of the system degrees of freedom x{s). This leads naturally to the "doubling" of the degrees of freedom 
in a real time description. 

Typically one assumes that the initial density matrix pxuQ{ti) is uncorrelated in its variables x and q, i.e. 

PxuqiU) = px{xt,x'^,U) (g) pQ{q^,q[\U), (A5) 

which can be motivated by assuming that the interaction >Cint(a;, q) is adiabatically switched on at the time ti. 
Introducing the reduced density matrix pr as 

Pr{x, x; t) = Tr(^) {p{x, q; x' , q; t)} = px{x, x'; t), (A6) 

one finds that the expectation value of any operator A depending solely on the system degrees of freedom can be 
readily expressed in terms of pr: 

(A) = Tr(^^^) {A{x, x'; t)pxvQ{x, q; x' , q'; t)} = Tr(^) {A{x, x'; t)pr{x, x'; t)} . (A7) 

Moreover one finds that its evolution in time can be put in the general form 

Prix,x']t) f dx,dx', f Vx f 2?x' e^(^=^[^l-^-[^'l' e^^'^[^'^'Vr(a;»,a;^;ii) (A8) 



by introducing the influence functional S'iF[a;,a;'] [2^ 

giS:p[..x'] ^ T^^^^ {u^Q^\q,q,;t,U)pQ{q,,q',-U)U^Q''^\q',,qf,t,U)} 

dqdq.dq'i P Vq f\q' e^^^'3l'l]-SQlq']+S,ntlx,q]-S.ntlx\q']) ^^f^^^^ ^/^.^^-^^ (^g) 



Here x{s),x'{s) are treated as given, classical background fields. Thus one is led to say that Pr evolves in time 
according to the effective interaction 

S,s[x,x'] = S.Ax] - S.,[x'] + Sif[x,x']. (AlO) 

It is clear that Pr{t) evolves causally from the history of the system and the bath. 

A notational simplification is achieved by introducing combined variables Xc{s), qds) defined on the real time 
contour C in the real time Green's function approach to finite temperature quantum field theory or nonequilibrium 
quantum field theory, i.e. 
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/ N _ J x{s) s £ upper branch / att\ 

'^^ ' 1 x'{s) s e lower branch, ^ ' 

and similarly for q{s), illustrated in Figure Al. The influence then takes the more compact form 

J J qi J q'. 

where the integration is defined on the time contour from ti forward to t and back to ti. 
For the influence action 5if[2;, x'] one finds the following general properties to hold: 

Sif[x, x'] = - {Sif[x', x])* , (A13) 

Sif[x,x]=0. (A14) 

Expanding Sif up to second order in x and x' , the general structure of Sif is given by 

Sif[x,x']k J ds F{s){x{s) - x'{s)) 
1 /■* 

+ - / dsids2{x{si) ~ x'{si)) R{si,S2) {x{s2) + x'{s2)) 

i /■* 

+ -J dsids2ix{si)-x'{si))I{si,S2){x{s2)-x'{s2)). (Al5) 

Suppose that Ci^t takes the quite general form 

Cint[xc,qc] = / {xc{s)) S (gc(s)) • (A16) 
Writing for Sif [x] up to second order in / 

exp —Sif[x ] « [•^'^ dsiF{si)f{x^{si)) + ^ dsids2fix^{si))D{si,S2)f{xa{s2))] (A17) 



it follows by comparing with (A12) and performing the substitution x{s) f{x{s)) 



i) 5f{xc{si)) 



^(^1) = ( 7 ) UTTTT^ - (S(<7c(si))),« 



D{si, S2) 



i J 6f{xcisi))Sf{xi{s2)) 



= [(P(S(g,(.i))Sfe(s2)))> - (Sfe(.i)))(S(gc(s2)))] (A18) 

where P means path-ordering along the contour C. F{s) can be interpreted as an external force term due to the mean 
field generated by the average interaction of the bath variables Q with the system X. 

The path-ordering definition leads to the typical four Green's functions defined in real time: 

D++{t,M) = \m^{q{hmq{t2))) - {^{q{h))){^{q{t2)))) 

= e{h - t2)D-+{h,t2) + e{t2 - h)D+-{hM) 

D+-ih,t2) = ^mqit2mqih))) - {Eiqih))){Eiqit2m ^ D<{t,,t2) 
D-+{t,,t2)^\mq{timq{t2))) (S(g(ii))>(S(g(t2)))) ^ D>{t,,t2) 
(il,i2) = J((f(S(g(ti))S(9(t2))> - {^{q{h))){E{q{t2)))] 



D 



= 6{ti ~ t2)D+-{h,t2) + e{t2 - h)D-+{t2,t2) (A19) 
Here T (T) stands for (anti)-timc ordering. 
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Defining the real Green's functions 

1 



Di{tiM) - - {D>{h,t2) + D<{ti,t2)) 



(A20) 



where Dfi{ti,t2) — —Dji{t2,ti) and Di(ti,t2) — Dj(t2,ti), one finds after some algebra that in ( A15| ) i? and / are 
given by 



R{ti,t2) = DR{ti,t2)0{ti-t2) = 2Re{D++{ti,t2))6{h-t2) 
lihM) = ]^Di{ti,t2) =lni{D++{h,t2))- 



(A21) 



If the system behaves quasi-classically, the density matrix pr becomes nearly diagonal. Then from (AS), the major 
paths contributing to the evolution of the reduced density matrix are obtained by extremizing the effective action 
Scs[x,x']: 



and the extremity conditions 



x{s) = x' {s) for all ti < s <t 



Sx{s) 



= 0. 



(A22) 



(A23) 



The equations of motion are clearly causal if one constructs the effective action to the argument s = t. 
Typically one introduces the variables 



X + x' 



A = 



so that the semiclassical equations of motion (A23) for x are stated as 

SScsix,A) 



SA{s) 



A=0 



From (A15) we have by expansion up to second order 

'SSx[x] 



thus 



n n 



SSx\x] 



Sx{s) 



Sx{s) 



F{s) 



F{s) + Ri^ X ] - —A (8) /(g) A, 



ds'R{s,s')x{s') = 



(A24) 



(A25) 



(A26) 



(A27) 



as the average quasi-classical equation of motion. This equation, in return, has to be interpreted as an average 
over random, fiuctuating forces. The last contribution in ( A26| ) leads to decoherence because any path contributing 
with sizeable |A| > over past time becomes exponentially suppressed — hence, the imaginary part / drives the 
system to quasi-classical behavior. For short periods in time, however, fluctuations in A can appear on the order of 
I / dtidt2l{ti,t2)\~^^'^ stochastically. These act as random "kicks" on the actual trajectory and can be interpreted as 
a stochastic force llTtl. To see this in more detail, one defines the real stochastic influence action 



SiF[x,x',£,]=Re{SiF[x,x']) + J ds ^{s) {x{s) ~ x' (s)) , (A28) 
where ^(s) is interpreted as an external force, randomly distributed by a Gaussian distribution with zero mean: 

(A29) 



^K(s)] = ^exp (~^J^ ds^ds2^{s,)I-\s,,S2)^{s2) 
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The influence functional S'iF[a;, x'] is regained as the characteristic functional over the average of the random forces. 

e'^'^^^^''"^ = (l)£_[x{s) - x'{s)] = I Di P[S,{s)] e*SiF[a;,x',C]^ (^30) 



The imaginary part / also defines the correlation of the random forces 

(C(siK(s2))5-?i/(si,S2). 



From (A28) the corresponding equation of motion reads 

5Sx[x 



5x{s) 



F{s) 



ds R{s, s')x{s') = ^(s) 



(A31) 



(A32) 



as the effective Langevin like equation dictating the dynamics of the system variable in the quasi-classical regime. 
(Classical) Brownian Motion jl^] is recovered if the spectral function of the system allows for approximating the 
kernels R and / as 



R{s, s') = —m'~i5{s, s') , 
hl{s,s') — 2ni'ykT5{s, s) 



(A33) 



In summary, by coarse-graining over th e bath d egrees of freedom, there appear two major contributions to the 
evolution of the reduced density matrix ( A8 JA26 ), namely the oscillatory or quantum mechanical phase and the 
decoherence factor when summing over paths in the (a; — A)-plane. We do expect an essentially semiclassical evolution 
when the excitation energy of the system (divided by Ti) is large so that the oscillatory phases of paths neighboring 
the classical path interferes destructively. Moreover, all paths where the time averaged A is sufficiently large will be 
exponentially suppressed by the decoherence factor. Hence, the system will be driven by decoherence into the quasi- 
classical regime, depending crucially on the magnitude yfi^^ ~ ^/l of the random forces, i.e. the size of the noise 
kernel (A31). If it is large, decoherence should happen very rapidly, and the quasi-classical description is inherently 
stochastic. 



APPENDIX B: THE INFLUENCE ACTION FOR SCALAR (A^ THEORY 



Dividing the scalar field <& into soft m odes </> and hard modes ip, the completely analogous expression for the 
influence functional can be read off (A12): 



J J {pl J ip'. 



where in the definition of tp] 



also the self interaction of the hard modes is included and shall be treated perturbatively. 
Expanding the influence functional up to order g'^, 

c _ o(0) I , 0(2) 

OIF — + + •JjF I 

one finds in an obvious notation: 



0(0) _ -, 



^(2) 

'if 



1 i 

2h 



(Bl) 



(B2) 



(B3) 



(B4) 
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The expectation values can be evaluated perturbatively with standard real time Green's function techniques p2| , p3[ . 
To avoid difficulties with initial higher order correlations of ph{ti) one typically has to assume that phiU) is of the 
form 



PQ;0 = Ph{U) exp Apfltop^ , {Ap < 0) 



(B5) 



to allow for a simple Wick decomposition. A thermal distribution of non-interacting hard particles is thus appropriate. 
We define the real time propagator of the hard particles on the contour as 

G','{xi,X2)^t{P{^{xiMx2)))o (B6) 

with its components in real time 

G'^^^{xi,X2) = i{T{ip{xi)(p{x2)))o = 9{ti ~ t2)G'!L%ixi,X2) + 0{t2 ~ ti)G'^^_{xi, X2) 
G'1^_{xi,x2) = i{(p{x2)yyixi))o = G'^{xi,x2) 
Gt%{xi,X2) = i(<^(a;i)<^(a;2))o = G'^^{xi,X2) 

G'lUxi,X2) = i{fiipixi)ip{x2)))o, (B7) 
which have a simple Fourier representation 

(p, r) ^ G^<- (p, -r) = ^ + '"("^^^ " ' for b| > fc. 

2ujp 

G^=(p,t) =G^ip,-r) =0 forH<fc, (B8) 



The evaluation of the expectation values in (B4) is lengthy but straightforward. We state only the final results with 
their real and imaginary parts: 

2 rC 

5}p"^=i^/ d^xG^,-ix,x)(bix)^ 

= ^^l' d'x im - 4>'{x)) G^>^ (0) {4>{x) + </>'(x)) (B9) 

^If = d^xrd^X2cj}{x^f (G^(xi,a;2))'G^(a;2,X2) 

d'*'Xid'^X2 {4>{xi) - (t>'{xi)) 

ti 

X \e{ti - t2) (Gl^{xi,x2f - G^^(xi,a:2)') Gi^iQ)] (<^(a;i) + <^'(xi)) (BIO) 

Sg^^ = d''xid^X2(t>{xi) (G^(X1,X2))'0(X2) 

= -Y^ / d^xid^X2{{(ty{xi)~cj,'{xi))e{h~t2) 



ti 



X (g^=(xi,X2)=^-GJ;=(xi,X2)3) (0(x2)+0'(x2))} 



X ^ (G^=(xi,a:2)'+Gi;=(xi,X2)3) (0(X2)-0'(X2))| (Bll) 



S[l''^ - -^f^ d^x,d^X2cl>{xif (G^(a;i,X2))%(x2)2 
16 
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-^f^ f d^x^d^X2{{4>{xif-<t>'{x^f) 



X i(G^=(xi,X2)^ + G^=(a;i,a;2)2) (0(x2)2 - 0'(a;2)') 
d** a; 1 d"* X2 (/) (x 1 ) ^ G^" (x 1 , X2 ) (x2 ) 



(B12) 



= r '^'^i'^'2;2 {{<t>{x,f - c^{x^f)9(tr - t2) 



(G^=(xi,X2) ~ G^=(xi,X2)) ((/.(X2)3 + 0'(X2)3)} 



+,|! f d^xid^X2{{Hx,f-cb'{xif) 
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^ (g^=(xi,X2) + G^=(X1,X2)) ((/.(X2)3 - (t,'{x2f 



(B13) 



Momentum space expressions can be found by inserting (B7) into these equations. Introducing the new field 
variables 



(B14) 



and using the relations 

Re(G^V(xi,a;2)") = {e{h - t2) + (-l)"0(i2 - ti)) ^ (Gi=(xi, 0:2)" + (-l)"Gj;=(xi,a;2)") 
Im(Gt+(xi,a;2)") = - i2) - (-l)"0(t2 - ti)) ^ (g^^ xz)" - (-irG^;^ (xi, a;2)") 

one finds after some calculation for the real and imaginary part, respectively, of the influence action 



(B15) 



Re S'if[(?!'a,0] = — ^ J d'^xi(j)Aixi)(t>{xi) Im G+'=+(0) 



dVd*a;20(ti-t2)|2<^A(a;i)</.(xi)ImG;;^+(O)Im G^=+(a;i, xs) 
"■^0a(xi)Rc G'l\{xi,X2f (t>{x2) 

+ ^4>{xi)(I>a{xi) Im G'l\{xi,X2f (^^{X2)^ + ^(t)A(x2f 
+ l<pAixi) U'(xi) + ^0A(a^i) j Re [g';^+(xi,X2)] ^(xz) U(x2)' + ^(bA{x2f 



(B16) 



and 



Im S*] 



IF \(PA, 



d'*xid''x2 |-— 0a(xi) Im G'i%{xi,X2)^ <Pa{x2) 
-^0a(xi)0(xi) Re G5ff+(xi,X2)^ 0a(x2)0(x2) 
+ ^0a(xi) U(xi)^ + ^0a(xi)2 



xlm 



G'J_"+(xi,X2) (/)a(x2) ( 4>{x2)'^ + J^4'a{x2)'^ 



(B17) 



After introducing the stochastic auxiliary functions ^j(x), z = 1,2,3, the imaginary part of the influence action 
becomes 
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In S'if[0a,0] — > •S'lF [0, 0A;^i] 



The corresponding stochastic weights are 



P[ii] = Ni { / d'xid^X2£,i{xi) 



Im 



PM =N2\ d^Xld^ X2UX1) 



9 

1 r 



Rc(Gl 



12 



6(^2) 



6(2^2) 



Im g: 



UX2) 



(B18) 



(B19) 



APPENDIX C: DERIVATION OF TWO DAMPING COEFFICIENTS 

Following the interpretation developed in Section III. A, the dissipative contribution due to Re ^S'lp'^ to the equatic 
of motion of the low momentum mode (/'(k, t) is stated as: 



'^■^ d^k2 



(2^)3 T(l(k-k2,<)^(k2,i) 



d^h 



i2n) 



' T(l(k + k2,t)0*(k2,t) 



'—^ <^(k - ki + k2, i) 0(ki, t)</)*(k2, t) 



1 



X I iM^'''>{k + k2, £;k-ki+k. +Ek,)+ iM^'^Hk + k2, Sk-ki+k. - Ek, 



0(fc,-|k-ki+k2|) 



(CI) 



This is a highly nonlinear term where the integrand rapidly oscillates over the whole integration region. However, its 
'major' contribution is expected for momenta ki « k2 because then ^(ki, t)(/)*(k2, i) « |(/)(ki,i)p > 0. It is not our 
objective to show that all other contributions cancel, but to extract the dominant contribution to the integrals in (|C1|). 
Indeed, setting ki = k2 corresponds to the Wick contraction in standard thermal field theory. Suppose the above 
expression would be taken as an ensemble average over all possible and similar configurations, i.e. initial conditions, 
thus averaging over all the fluctuating phases. For a quasi-homogeneous system one will find for the correlator 



Here 



(^(ki, i)0*(k2, t)) « (27r)-^5^(k2 - ki) |</>(ki, i)| 



|0(ki,i)|2«-L7V(ki,t)« -L(l + iV(ki,t)) , 



Ek, 



(C2) 



(C3) 



which is identified with the average population number in the classical regime, when the mode population is large. 
(C2) is then completely analogous to the equal-time limit of the propagators (defined in Appendix B) in a homogeneous 



system: 



iG>(k, r = 0) « n{uJk)/uJk . 



(C4) 



As a further motivation let us first calculate the classical vertex contribution to the equation of motion by means of 
the approximation (C2): 



ff2 d^kid^k2 



3! 



g (/)(k-ki -k2,i)0(ki,i)0(k2,t)0(fcc- |k-ki -k2|) 

|0(ki,t)p ^(k,t) = Am2i0(k,t) . 



(2^) 



(27r)3 



(C5) 
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Neglecting the fluctuations let us write the average influence of the other soft modes as a mass term induced by the 
mean field. This mass term corresponds to the lowest order tadpole contribution in the quantum field theory. Noticing 
the identification (|C^), we find that Am^j exactly equals the expression (^) stated in the Introduction. Now, writing 

</>(k- ki +k2,i)'/'(ki,0'/'*(k2,t) = 0(k,t) {(t){ki,t)(j)*{k2,t)) + fluctuations , (C6) 



the first term contributes to ( |Cl[ ) as: 

d^ki 



'/'(k, i ) — y ^ I </)(ki , t) 1 2 (zX ('^) (k + ki , i?k + J + * ('^Hk + ki , iJk - £;fe J 

= 77W(k,t)0(k,i) . (C7) 



With the identification (C3) this expression thus corresponds naturally to the scattering diagram obtained by cutting 
the 'sunset' diagram with two hard and one soft propagator as illustrated in Figure CI. 

Similarly, for the dissipative contribution of Re ^S'jp''^ we write for the first term: 

^!M^^!||jl^^(k3, 0(ki, t)0(k2, i)0*(kl, t)</.*(k^, t) 0(fc, - |k3|) 

x-|^ (iM^''\k + k'l + k^, ^1 + ^2 + S3) + iA^("'(k + k'l + k^j, -Si - ^2 + £^3) 
0-C/3 V 

+227W(^)(k + k'l + k^, Si -E2 + S3)) (C8) 

with k3 = k + k'l + kj — ki — k2 . The other term is of higher order in the time derivatives and can be neglected if 
the mode population changes weakly in time, i.e. 

||0(k,i)p « Sfc|</.(k,i)|2 . 

Again, the dominant contribution to the integral arises when ki k, k'l, k2 ~ kj or ki « kj, k2 ~ k'l, so that we 
approximate for a quasi-homogeneous system: 

(0(ki,t)0(k2,i)r(k'i,t)0*(k^,t)) 

« (2^)6(^3(k^„k'i)^3(k2„k'2) + 53(ki-k'2)<53(k2-k'i)) |0(ki,t)nc/,(k2,t)P (C9) 

Inserting this into ( |C8[ ) one obtains 

<^(k,i) ^1'^ i<^(ki,i)n<^(k2,op 

X [iMS"^ (k + ki + k2, Si + S2 + Sfe) + iMS"^ (k + ki + k2, -Si - S2 + Sfe) 

+ 2zX('^) (k + ki + k2, Si - S2 + Sfe)) 

= 77(^)(k,i)(/)(k,t) . (CIO) 

This corresponds to cutting the 'sunset' diagram with two internal soft propagators and one hard propagator giving 
rise to the various on-shell scattering and production processes as illustrated in Figure C2. 
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(a) (b) (c) (d) (e) 

FIG.l. Feynman diagrams contributing to the influence action up to order g*. Full lines denote hard modes and dashed lines correspond 
to soft modes. 



Q... 



Re 



FIG. 2. Graphical representation of the classical equation ( ^ ) for the soft field modes. The noise terms 
with (N — 1) external legs. 
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(a) 



(b) 



FIG. 3. (a) A particular contribution to the effective action at order in <l?^-theory; (b) Lowest order diagram leading to on-shell 
dissipation when iterating the real part Re(5|'^') of (a) together with the ^'P'^ term in the full nonlinear and semiclassical equation of 
motion. It corresponds to the usual two loop diagram when cut through the hard thermal loop and the soft (thermal) field line. 
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x'(s) 

FIG.Al. Schwinger-Keldysh time contour path C for the variables x{s) and x'{s) running from s = tj to +00 and back to tj. 
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FIG. CI. Various scattering contributions when cutting the "sunset" diagramm with two hard particles and one soft mode. 
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FIG.C2. Various scattering contributions when cutting the 'sunset' diagram with one hard particle and two soft modes. 
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